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We analyze a 2D spin-pseudospin model, where the pseudospins represents the charge 
degrees of freedom. The model is known to undergo a phase transition with the simultaneous 
appearance of the long-range charge order and the spin gap. We show how the gap vs 
critical temperature ratio (also called the BCS ratio) gets renormalized from the classical 
non-interacting value. This value is also universal in the sense that is does not depend 
' on the microscopic parameters of the model, and must be the same for various types of the 

. Peierls-like transitions where the spin gap is accompanied by the structural, orbital or charge 

^ ' order. 

o ■ 

' • Introduction. ~ There is a large number of materials which demonstrate phase transitions into states 

with a spin gap accompanied by some types of the structural, charge or orbital order. The most canonical 
example is the spin-Peierls transition [l|, [3] ■ The interplay of charge, spin and orbital degrees of freedom is 
^ . known to produce some exotic phases in transition metal oxides [3, [j] . The dimerized Peierls states driven 
by the superstructures of the orbital order are reported in some spinels Q]- Our main motivation comes from 
the recent work on the spin-SAF transition in the quarter-filled ladder compound NaV205 where the 

■ Super-Anti-Ferroelectric (SAF) long-range charge order occurs together with the spin gap. Several other 
layered vanadate compounds demonstrate transitions when the spin gaps occur simultaneously with charge 



o 



o 

(N 



' ordering. In particular, the spin-SAF transition was recently reported in Zn(pyz)V40io 



If, in all seemingly different Peirels-like states the spin gap is induced by the dimerization due to 
O . structural, orbital or charge long-range order, then there should be some universal parameters unifying 
all such transitions. A good candidate is the BCS ratio we calculate here for the case of the spin-SAF 
transition. It does not depend on the model microscopic parameters and matches the value found earlier for 
the spin-Peierls transition. 

• Spin-Pseudospin Model and Analysis. - We analyze here a spin-pseudospin model which consists of 



the Ising pseudospins T coupled to the Heisenberg spins S which reside on the same sites of a square lattice. 
The pseudospin sector is given by the Ising Model in a Transverse Field (IMTF). An elementary plaquette 
and the couplings in this model are shown in Fig. [TJ The Hamiltonian of this IMTF reads: 

X. 

= The spin sector is given by the Heisenberg chains parallel to the Ji diagonals Hs = J J2m n ^mnSm,n+i- The 
spin-pseudospin coupling is: 

HsT = SmnSm,n.-|-l (7^+l,n+l — 7^-l,n) ' (2) 

m,n 

with the total Hamiltonian of the model: 

H = HiMTF + Hs + HsT ■ (3) 

The spin and pseudospin operators satisfy the standard SU (2) algebra, while S and T commute. 

The spiri-pseudospin model ([3|) and some of its modifications were proposed and analyzed in the earlier 
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related work [THlOl] in the context of the quarter- filled ladder compound NaV205, where the pseudospins 
correspond to the charge degrees of freedom. For the physically interesting couplings [lOjjthe classical Ising 
model {Q = 0) orders into the 4-fold degenerate Super- Anti-Ferromagnetic (SAF) phase shown in Fig. [TJ 
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FIG. 1: Couplings on an elementary plaquette of the 2D nn and nnn Ising model (left) and an example of the 4- 
fold degenerate SAF order (right). The SAF is also called columnar or stripe order in some recent literature. The 
degeneracy of the SAF state is due to the Z2 ® Z2 symmetry with respect to the spins of each of the sublattices (labeled 
by circles or squares) being flipped. The lattice vectors are parameterized by the integers (to, n) in the skewed basis 
r = TOBi + ne2. 



The SAF pattern appears in various contexts, and it is also called columnar or stripe order in some more 
recent literatm'e. The SAF charge order occurs along with dimerization and gap in the spin sector lOl ]. 
That is what we called the spin-SAF transition Note that the spin gap is due to the frozen phonon 

displacements at the spin-Peierls transition, while the charge plays the role of phonons at the spin-SAF 
transition. 

The molecular-field approximation is applied for the pseudospins [7,], while the spin sector is treated 
via minimization of the exact free energy of the dimerized Heisenberg XXX-chsin 



i^xxx = 5^j(i + (-ir5)s„s„+i 



(4) 



The specific free energy of the spin chain /s(T, 6) is an analytic function at T 7^ and can be expanded over 
5 as fs{T,5) = fs{T,0) — ^r]{T,0)6'^ + 0(5^), where r/(T, 0) is called the static dimerization susceptibility. 

• Spin-SAF phase transition. XY Spin Chain: Free Fermions. ~ It is straightforward to obtain in a 
closed form the specific free energy of the XY spin chain mapped onto the spinless non-interacting Jordan- 
Wigner fermions. To leading order 
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(5) 



where C = — ^r^, 7 = 0.5772... is Euler's constant. In the region J < Jc (where Jc = 2Q is the mean-field 



value of the QCP of the IMTF ([T]) [lO|), the critical temperature is given by the BCS-type solution 



Tc^CJexp -^(J^-J) 



The ground-state dimerization is 



6 w -exp 

e 
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(6) 



(7) 



In the case of free fermions the spin gap depends linearly on dimerization, A-'^'^ = J6. So the ratio of the 
zero-temperature spin gap (Aq) and the critical temperature (a.k.a the BCS ratio) in the regime J7c < J7 is 



A^ 
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e^ 



1.76... 



which coincides exactly with the classical result for the superconducting gap in the BCS theory [1 
• Spin-SAF phase transition. XXX Spin Chain: Interacting Fermions - sine- Gordon Model. 



(8) 



The 



Heisenberg spin chain can be mapped onto the model of interacting Jordan- Wigner spinless fermions, and 



the low-energy sector of the fermionic Hamiltonian in its turn can be bosonized 13(] . Neglecting the marginal 
term, the dimerized spin chain maps onto the sine-Gordon model 13| 



(9) 



where v = is the bosonic velocity and = ^6. The relevant perturbation of the free bosonic part of the 
Hamiltonian ([9]) comes from the spin dimerization term (— l)"'S„S„+i ~ cos \/27r(/). The amplitude is 
not known exactly yet, but according to the approximate calculations of Orignac 1G| 



(10) 



Using the sine-Gordon model ([9j) to approximate the low-energy sector of the dimerized Heisenberg chain 
dH) , the free energy of the latter reads to leading order [ij, [13] 



UT, 6) = -JU -i-- l^S' , where = ^ = ^.A^ 
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and to = ln2 — ^ = 0.4431.... Then the dimerization susceptibility to lowest order 

V{T,0) = ^. (12) 

This function was first calculated by Cross and Fisher ^ for their spin-Peierls transition theory. The 
amplitude A^ pH]) suggested by Orignac [la] gives ao = 0.253418, very close to the original bosonization 
result a^" ^ 0.2610]. 

The interactions of the Jordan- Wigner spinless fermions of the XXX chain modify considerably the 
critical properties of the Ising-XXX model The inverse dependence rj oc ()12p results in the power-law 
dependence of Tc on couplings: 



Jc-J 

The phase diagram of the Ising-XXX model is shown in Fig. [5J 



J < Jc 



(13) 




FIG. 2: Phase diagram of the coupled Ising-XXX model. Axes are not in scale. The solid red line shows the critical 
temperature Tc separating the disordered and SAF phases of the decoupled IMTF. Its QCP corresponds to Jc = 2f2 in 
the mean- field approximation. (J = Ji + J2). The dashed blue line shows Tc separating the disordered -f spin gapless 
phase from the spin-SAF phase of the coupled Ising-XXX model. 



From the ground-state specific energy of the sine-Gordon model 15|, [18] 



. X M2 vr , 2 / F 3/4 \^/^F l/6 

/sG 0, a) = -V tan - , where M = — ttu ) ' ( ' (14) 
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(M is the dhnensionless soHton mass) we get the ground-state energy of the Heisenberg cham 

/,(0,Ao) = -Jio-^^ , (15) 
where the zero-temperature spin gap Aq = ■^JM is related to the dimerization as follows: 

The ground-state spin gap: 

A,=^,A?^, where *..?,/^fn^yf£iVfiy (iv) 

'Jc-J 3\3\r{i/A)J Vr(2/3); ^ ' 

Combining Eqs. ()13ll7p we obtain the BCS ratio: 

^ = ^i = 6v/3iniM!.2,47.,. (18) 

Tc fi {r(i/4))' 

The same BCS ratio for the spin-Peierls transition was first obtained by Orignac and Chitra [l^. We would 
like to stress that the above result is exact for the IMTF coupled to the sine-Gordon model, i.e., when 
the marginal terms in the Heisenberg spin Hamiltonian are neglected. Similar to the non-interacting result 
dS]), the ratio (llSp does not depend on the microscopic parameters of the model, and even the dimerization 
amplitude cancels. In this sense we interpret this as a universal result. The interactions renormalize the 
BCS ratio away from the free fermionic value of 1.76. 

• Conclusions. -The BCS ratio in the interacting Ising-XXX model is calculated. Similar to the classical 
free-fermionic case, this value is also universal in the sense that is does not depend on the microscopic 
parameters of the model. We conjecture that it must be the same for various types of the Peierls-like 
transitions where the spin gap is accompanied by the structural, orbital or charge order. An extension of 
these results taking into account marginal terms of the spin chain Hamiltonian is warranted. 
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